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Abstract 

We consider the possibility of a direct determination of the polarized gluon 
density in the proton using charm production with polarized beams at HERA. 
We study total cross sections and distributions at leading order using different 
parametrizations of the polarized gluon density. We conclude that charm pho- 
toproduction data can be used to constrain the polarized gluon density in the 
proton if an integrated luminosity of 100 pb~ x will be achieved at HERA. 
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A direct measure of the spin-dependent gluon density in the proton, Ag(x,Q 2 ), 
has never been performed (see ref. |1] for a recent review on polarized nucleon struc- 
ture). Polarized deep inelastic scattering data allow to extract the structure function 
gi(x,Q 2 ) at different values of x and Q 2 and to obtain an indirect determination of 
Ag through scaling violation 0-|5j . However, such indirect determinations are affected 
by large uncertainties, because of limited statistics and limited coverage of the x, Q 2 
range of the presently available experimental data. On the other hand, a precise 
knowledge of the polarized gluon density is important in order to understand various 
aspects of the structure of polarized nucleons. For example, it has long been known || 
[7| that the anomalous gluon contribution to the first moment of g±, although formally 
of order a s , does not vanish in the large-Q 2 limit, and therefore affects the extraction 
of the singlet axial charge from polarized deep inelastic scattering data. Further- 
more, an independent determination of Ag would allow to test the reliability of the 
perturbative expansion, and to assess the importance of possible non-perturbative 
contributions. 

In order to measure Ag directly, it is necessary to consider processes which are 
only or predominantly initiated by gluons; many have been studied in the past, like 
for example heavy quarkonia production in photon-gluon fusion || , or the production 
of large-fey jet pairs in deep-inelastic scattering PUlQfl. Another interesting possibility 
is the production of heavy quarks in photon-proton collisions. This process has al- 
ready been considered in refs. [p7T] - p73j| in view of the possibilities of the electron-proton 
collider HERA, and the conclusion was reached that the study of total cross section 
asymmetries for heavy-quark photoproduction at HERA can be of little help in con- 
straining the polarized gluon density in the proton. In this letter we reconsider this 
problem, in the light of recent experimental information on polarized deep inelastic 
scattering (|14| , |T5|| and improved theoretical understanding, and of the planned values 
of luminosity for the HERA collider. We will also consider the impact of realistic 
experimental cuts on the observed quantities, and discuss some theoretical uncertain- 
ties. We will always refer to the production of charm quark-antiquark pairs, since 
this is the most favourable case for the kind of study we are considering here. 

The electron-proton cross section for the production of heavy quarks can be reli- 
ably computed in the Weizsacker- Williams approximation |TjJ, where only the con- 
tribution of on-shell collinear photons radiated by the incoming electron beam is 
retained. Alternatively, the HERA experiments are capable to tag the scattered elec- 
tron, thus determining the energy of the photon that initiated the reaction under 
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study. In this kind of analysis, HERA works as a photon-proton collider with a well 
defined photon energy, ranging between 70 and 270 GeV approximately In both 
cases, the relevant parton subprocess at the leading order is 

g(P,X 9 ) + i(qA 1 ) ^c(k) + c(k'), (l) 

where four momenta and helicities are indicated in brackets. The 0(a em a s ) partonic 
cross section for the process ([!]) is given by p"T| . |T2" 



6 GCq 

da 19 {s, t, X g , Xy) = -^r [S + A 9 A 7 A] f3dcos9, (2) 



where 



8m 4 s 2 i 2 + u 2 + Am 2 s 

+ ^ j (3) 



t 2 v? tii 

4m 2 (t 3 + u 3 ) i 2 + u 2 — 2m 2 s 

4 + 2 — 

t 2 u 2 tii 



*= l7 2 "' +2 ' ^ (4) 



m is the charm quark mass, e c its electric charge, (3 = yl -4m 2 /s and 6 1 is the 
scattering angle in the 7$ center-of-mass frame. We have defined 

s = (p + q) 2 , t — (p — k) 2 — m 2 , u = (p — k') 2 — m 2 . (5) 

The quantity of interest is the difference 

dAa w = ~ (da\l - da\ l p ) , (6) 

where dcdj and dcr^ are the differential cross sections for the cc photoproduction 
process, with parallel and antiparallel polarizations of the incoming photon and pro- 
ton, respectively. In fact, it can be easily shown that at leading order the quantity in 
eq. (|D can be written as 

dAa lp (s, t) = Ag(x g , fi 2 F )dAa ig (s, i, fi 2 R )dx g , (7) 

where 

dAa 7g (s, t) = - (da ig (s, i, +1, +1) + da 19 (s,i, — 1, -1) 

-da ig (s, i, +1, -1) - da 7g (s, i, -1, +1)) . (8) 
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In eq. (0) we defined s = (p + q) 2 and t = (p — k) 2 — m 2 , where p is the four-momentum 
of the incoming proton, and x g is the fraction of the proton longitudinal momentum 
carried by the gluon; fi R and fi F are the renormalization and factorization scales. 

In the case of electroproduction, using the Weizsacker- Williams approximation we 
have 

dAa ep (s, t) = Af ww (xry, Q 2 ww )Ag(x g , fi 2 F )dAa ig (s, i, ffydxgdx^ , (9) 

where 

Af ww (x„Q 2 ww ) = ^ ~ (1 ~ ^ log ^]-^\ (io) 



e 7 



and we defined s = (p+q) 2 , q being the electron four- momentum. The mass scale Q 2 WW 
entering the Weizsacker- Williams function has been chosen as discussed in ref. |T7| ]. 

We will consider the three fits to Ag presented in ref. p|, which we will indicate 
with GS-A, GS-B, GS-C, and those presented in ref. Q, denoted by BFR-AB, BFR- 
OS and BFR-AR. The three parametrizations of Ag given in ref. [|J are obtained 
by performing fits to data within three different subtraction schemes for collinear 
divergences. Since a next-to-leading order calculation of the polarized partonic cross 
sections is not available, the scheme choice is immaterial in our analysis. 

We begin by considering total cross sections. In the first row of table [I] we present 
(columns indicated with I) the total charm production cross section computed at 
next-to- leading order |l8| - |20l . We show both the electroproduction results, obtained 
using the Weizsacker- Williams approximation at the HERA center-of-mass energy 
of 314 GeV, and the photoproduction results, with a center-of-mass energy of the 
photon-proton system of 200 GeV. The factorization and renormalization scales are 
chosen equal to 2m and m respectively (see ref. [21] for a detailed discussion of scale 
choices) . 

We also show (columns II) the values of total cross sections with the conditions 

p T > 2 GeV (11) 

M < 1-5, (12) 

imposed on the transverse momentum p T and the pseudorapidity rj of the observed 
heavy quark, (assuming that only one of the two heavy quarks produced is fully 
reconstructed). The conditions ( ITT] ) and ( |T2"D approximately reproduce the present 



experimental situation of the HERA experiments. Finally, we show (columns III) the 
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ep 


IP 




I 


II 


III 


I 


II 


III 


a (nlo) (fib) 


0.5505 


0.09653 


0.04542 


3.929 


0.9205 


0.5268 


Act/ct gs-a 


7.19- 10~ 4 


5.66 • 10~ 3 


9.58 • 10~ 3 


2.78 • 10~ 3 


1.87- 10- 2 


2.82- 10~ 2 


Act/ct gs-b 


5.97- 10" 4 


5.96 • 10~ 3 


9.80- 10~ 3 


2.51 • 10" 3 


2.04- 10~ 2 


2.98- 10~ 2 


Act/ct gs-c 


3.87- 10" 5 


3.56 • lO- 3 


5.41 • 10~ 3 


3.38- 10" 3 


1.34- 10~ 2 


1.79- 10- 2 


Act/ct bfr-ab 


5.01 ■ 10" 4 


8.04 • 10~ 3 


1.28- 10~ 2 


2.03 • 10" 3 


2.87- 10~ 2 


4.02- 10" 2 


Act/ct bfr-os 


5.71 • 10" 4 


5.56 • 10~ 3 


9.20- 10~ 3 


2.24- 10" 3 


1.91 • 10~ 2 


2.79- 10~ 2 


Act/ct bfr-ar 


5.49- 10~ 4 


6.46 • 10~ 3 


1.05- 10~ 2 


2.26- 10" 3 


2.26- 10~ 2 


3.24- 10" 2 


1/V2a£ 


9.53- 10~ 5 


2.28 • 10~ 4 


3.32- 10~ 4 


2.77- 10~ 4 


5.73- 10~ 4 


7.57- 10~ 4 



Table 1: Total cross sections and total cross section asymmetries for cc 
production in ep collisions at ^fs^ = 314 GeV and in jp collisions at ^/s lp = 
200 GeV, for different choices of the polarized gluon density. The charm 
quark mass is 1.5 GeV, and the unpolarized gluon distribution is MRS A. The 
integrated luminosity for the ep system is 100 pb^ 1 . 



effect of applying the Peterson fragmentation function [^2J to the produced charm 
quarks, in order to estimate the impact of hadronization phenomena. In columns III 



the cuts of eqs. (|TT|) and (|12D are still applied. Notice that in the absence of the cuts 
the total cross section values would be insensitive to fragmentation effects. On the 
other hand, the cuts introduce in the total cross section a dependence upon the p T and 
rj distributions; in particular, since Peterson fragmentation softens the p T spectrum, 
the low-p T region gives a contribution to the cross section larger than in the bare 
quark case. This explains the lower values of columns III with respect to columns II. 

In the following six rows we display the values of the asymmetry Aa/a in the same 
cases, obtained with the six different parametrizations of Ag. It must be stressed that 
the unpolarized cross section a that appears in the denominator of the asymmetry 
is computed at the leading order, because a next-to-leading order calculation for the 
polarized cross section is not available. However, one might expect that the effect 
of radiative corrections approximately cancels in the ratio Aa/a. Notice that the 
asymmetries obtained with the fragmented quarks are larger than those for bare 
quarks, at variance with the case of the total cross sections. This is due to the fact 
that the dominant contribution to Aa/a is given by the large-p T region, where the 
asymmetry values are larger for fragmented quarks than for bare ones. We will discuss 
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again this issue when dealing with differential cross sections. 

The next-to-leading order value of a can then be used to estimate the sensitivity of 
the experiment. A rough estimate of the minimum value of the asymmetry observable 
at HERA can be obtained by requiring the difference between the numbers of events 
with parallel and antiparallel polarizations of the initial state particles to be larger 
than the statistical error on the total number of observed events. This gives 

"Act" 

(13) 



a 



n 

where £ is the integrated luminosity and the factor e accounts for the experimental 
efficiency for charm identification and the fact that the initial beams are not com- 
pletely polarized. The values of 1/ \j2o£ in the various cases are given in the last row 
of table [I], for £ = 100 pb _1 . In the case of photoproduction, the effective luminosity 
has been estimated by assuming that the results obtained with photons such that 
170 GeV< ^fs^, < 230 GeV can be reliably described by the cross section for an in- 
coming monochromatic photon with ^/s lp = 200 GeV. The ep luminosity of 100 pb _1 
has then been rescaled by a factor equal to the integral of the Weizsacker- Williams 
function in the appropriate range. 

By inspection of table [l], we conclude that the asymmetries of practical interest 
(columns II and III) are always larger than the corresponding minimum observable 
values with an integrated luminosity of 100 pb _1 even if an experimental efficiency 
of 1% is assumed. We also observe that the asymmetry is generally larger in the 
phase space region defined by eqs. ( |TTD and (|12]). In practice, the quantity £Aa is 
required to be larger than 5 or 10 events, in order to distinguish the signal from the 
background. This condition appears to be fulfilled for all the cases considered in 
table p], with £ = 100" 1 . We also observe that all the chosen parametrizations for the 
polarized gluon density give similar results for the total cross section asymmetry. It 
seems therefore difficult to distinguish among them with this kind of measurement. 

We now turn to differential distributions. From table [I] we know that the asym- 
metry tends to become larger in the phase-space region where the conditions (|ll]) and 
(12) are satisfied; we will therefore present our predictions with the kinematical cuts 



IT]) and (y^) applied. We will consider only the GS-A and BFR-AB parametrizations. 

In fig. || we show the asymmetry versus the transverse momentum of the observed 
quark for ep collision at a center-of-mass energy of 314 GeV. The scale choice is 
— /^O; fJ'F = 2/i , with /x — \/Pt + m2 - As previously observed, the asymmetry 
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increases with increasing transverse momentum. Observe also that the asymmetry for 
fragmented quarks is always larger than that for bare quarks over the whole p T range 
considered, consistently with the analogous behaviour of total cross sections. The 
results obtained with the two parametrizations for Ag are quite close to each other, 
although they display slightly different shapes: the BFR parametrization produces 
larger asymmetries in the \ow-p T region. The effect of Peterson fragmentation on 
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Figure 1: Asymmetry cross section versus transverse momentum in ep col- 
lisions at = 314 GeV. The minimum observable asymmetry, computed at 
next-to-leading order, is also displayed. 

the shape of this distribution is very small. In fig. p] we also show the quantity 
defined in eq. ([13]) for e = 1 and £ = 100 pb _1 , as an indication for the minimum 
observable value of the asymmetry. In this case, the quantity a appearing in the 
RHS of eq. ( |I3"D is the next-to-leading order cross section for a given p T bin. We 
used a bin size of 1 GeV; clearly, a larger bin size would decrease the minimum 
observable asymmetry; on the other hand, by enlarging the bins the resolution of 
the measurement would get worse. From fig. [I], we can see that there exists a wide 
p T region where the predicted asymmetries are more than one order of magnitude 
above the corresponding minimum observable asymmetry. Furthermore, the value of 
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the asymmetry for p T above 10 GeV is of the order of few tens of percent. We can 




x K, O; l/V 2fLo-(NLO) : 

H L=100 pb" 1 , f=0.0165 " 

ok y 



■_ K K; bare, <>: fragm. 

i i I i a i i I i i i i I i i i i 
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Figure 2: Asymmetry cross section versus transverse momentum in jp col- 
lisions at = 200 GeV. The minimum observable asymmetry, computed at 
next-to-leading order, is also displayed. 

therefore conclude that the p T distribution could be used to extract information on 
the polarized gluon density, even with a low experimental efficiency The shape of the 
distribution could in principle provide with the possibility of discriminating among 
the various sets for Ag. However, since the predictions for different sets are quite 
close to each other, this kind of measurement does not appear to be feasible with the 
planned luminosity values. 

Figure is the analogous of fig. [I] for 7p collisions at a center-of-mass energy of 200 
GeV. The asymmetry values turn out to be larger than for ep collisions in the whole 
p T range, and the curves appear to be slightly flatter than in the previous case. On 
the other hand, the minimum observable value is larger than before. We computed 
the effective luminosity for the incoming photon as in the case of total cross sections. 

It is also interesting to consider fully exclusive distributions, despite the fact that a 
large statistics is needed in order to reconstruct completely both the produced charm 
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and anti-charm. In analogy with the unpolarized case (see ref. ||23||), we can rewrite 
eq. (H) as 

' ,A<T,/ ' " -Af ww (x^Q 2 ww )Ag(x g ,fi 2 F )Aa yg (M^, (14) 



with 



dy c5 dM% s 



x 7 = — ™exp(-y C c), (15) 
V s 

x g = -y^exp(?/ c5 ). (16) 
V s 

In eqs. (0), (|15D and ( |16|) M^ E = s is the invariant mass of the cc pair, and y c5 its 
rapidity in the electron-proton center-of-mass frame. We have positive rapidities in 




Figure 3: Asymmetry cvoss section versus Xg in ep collisions at ^/s e p — 
314 GeV. The minimum observable asymmetry, computed at next-to-leading 
order, is also displayed. 



the proton direction. Identifying the LHS of eq. ( [14] ) with the experimental data, we 
can invert this equation to get 

(j . \ data 

dy c5 dMi 5 ) ■ { ) 



-9- 



We present our results directly in terms of the variable x g , exploiting the identity 



dAa ep dAa ep 

Xq 



dy c5 dM* 5 9 dx g dM%' v ; 

In fig. § we show our predictions for the x g asymmetry in ep collisions at a center-of- 



mass energy of 314 GeV. In this case, the reference scale /j,q is equal to y (p1 t + P§ T )/2 + m 2 . 
We do not present the result for the fragmented quarks. In fact, when Peterson frag- 
mentation is applied, x g is not expressed in terms of the invariant mass and rapidity 
of the pair as in eq. fll6|) . The correct relationship could in principle be worked out, 
but it is more straightforward to obtain the cc cross section by directly performing the 
deconvolution on the measured cross section for charmed hadron production. This 
procedure can be applied using different hadronization models, thus obtaining an es- 
timate of the dependence of the result upon the assumed hadronization mechanism. 

By inspection of fig. |3], we can conclude that the polarized gluon density can 
be probed with a good resolution (we used a bin size of 0.2 in the variable logx g ) 
in the range (10~ 2 , 10 _1 ), if the experimental efficiency is not too small. The GS- 
A and BFR-AB results are quite similar, and it appears unlikely that this kind of 
measurement could discriminate between them. 

In this paper we have studied the problem of measuring the polarized gluon density 
in the proton using the charm data which will be collected at the HERA collider in the 
polarized configuration. We used a leading-order QCD calculation, since the polarized 
short distance partonic cross sections are not available at the next-to-leading order. 
For this reason, we presented our results in terms of asymmetries, which are more 
likely to be stable under radiative corrections than the polarized or unpolarized cross 
sections. We have shown that, assuming an integrated luminosity of 100 pb _1 , the 
total cross section asymmetries can be measured even with a very low experimental 
efficiency. The measurement of single-inclusive and double-differential distributions 
also appears to be feasible, although in the latter case a higher luminosity is certainly 
needed. With a luminosity of 100 pb _1 , it seems difficult to distinguish among the 
various parametrizations for the polarized gluon density. 

We have not considered the inclusion of the hadronic component of the photon, 
which is known to be important in the HERA regime. Unfortunately, no information 
on the polarized parton densities in the photon is available up to now, and therefore 



no theoretical estimate can be given. It is known [|24] that the p T spectrum of the 
hadronic component is softer than the pointlike one in the unpolarized case, and 



therefore the cut of eq. (O]) reduces the impact of the hadronic component on the 
full unpolarized results. This might not happen in the polarized case. However, for 
hadronic photons a sizeable fraction of the photon momentum is lost into hadronic 
fragments, that can be observed in the photon direction. In this way it should be 
possible to disentangle the hadronic component from the pointlike one. 
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